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PROCEDE D'EXAMEN RADIOLOGIQUB D'UN OB JET, COMPRBNAWT TOl 
ETALONNAOE DB S PROPRIgTES D ' ACT - giroATION DTJ R AYONNEMBN T 



DESCRIPTION 



Cette invention concerne un proc&d£ 

d'examen radiologique d'un objet, oompronant 

e fealonnago deo preprint 6e d * as - t 6 nuafr ion c U* raLyonnemenU 
Lea proc^d^s radiologiques consistent & 

10 faire traverser un objet §. etudier par un rayonnement 
pour en d£duire la repartition de diffSrentes 
categories de materiaux, absorbant le irayoimeTnent ..de 
fa^ons dif£6rentes, dans cet objet, Une application 
tres courante est 1 1 oat 6odensitotnetri e , ou on recherche 

15 la massa ct la- density des tissus osseux dans un 
patient en distinguant ces tissus des tissus mous . 

II est usuel d'utiliser un spectre ;cle 
rayonnement large et & le diviser en bandes mesurSee 
s^paremant par des canaux de mesure respect if s, Cornme 

2 0 les coefficients d f absorption ou d 1 attenuation du 
rayonnement par l'une ou 1 • autre des categories de 
tissu sent differents pour chacune des bandes, le 
probleme thSorique se ramene £ la solution d'un systdme 
d' equations dont le nombre est 6gal £ celui. des bandes 

25 de mesure et qui comprennent chacune deux: iaconnues 



(les epaisseurs ou les masses traversees de tissus et - ^ 
osseux) . La solution du probldme devient possible aprds 
une calibration obtenue en faisant traverser par le 
rayonnement divers Stalons doc^a d'opaisseurs connues 
3 0 de matSriaux aux propriet6s d' absorption analogues B. 
celles des matferiaux de mesure, et notamment le 
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PROCEDE D'EXAMEN RAD I OLOG I QUE D 1 UN OB JET 
DESCRIPTION 

Cette invention concerne un procede 
5 d'examen radiologique multi-energie d'un objet. 

Les procedes radiologiques consistent a 
faire traverser un objet a etudier par un rayonnement 
pour en deduire la repartition de differentes 
categories de materiaux, absorbant le rayonnement de 
10 fagons differentes, dans cet objet. Une application 
tres courante est l'osteodensitometrie, ou on recherche 
la masse et la densite des tissus osseux dans un 
patient en distinguant ces tissus des tissus mous . 

II est usuel d ! utiliser un spectre de 
15 rayonnement large et a le diviser en bandes mesurees 
separement par des canaux de mesure respectifs. Comme 
les coefficients d 1 absorption ou d ' attenuation du 
rayonnement par 1 1 une ou 1 ' autre des categories de 
tissu sont differents pour chacune des bandes, le 
20 probleme theorique se ramene a la solution d'un systeme 
d' equations dont le nombre est egal a celui des bandes 
de mesure et qui comprennent chacune deux inconnues 
(les epaisseurs ou les masses traversees de tissus mous 
et osseux) . La solution du probleme devient possible 
25 apres une calibration obtenue en faisant traverser par 
le ^ rayonnement divers etalons dotes d' epaisseurs 
connues de materiaux aux proprietes d 1 absorption 
analogues a celles des materiaux de mesure, et 
notamment le plexiglas et 1 ' hydroxy apatite pour simuler 
30' les tissus mous et les tissus osseux. On peut alors , 
calculer les parametres mathemat iques d'un modele 
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plexiglas ec 1 'hydroxyapatite pour simuler les tiseus 

mous et les tissus osseux. L 1 inv ^ rcion d ' un oyoteme - 

d 1 equation — an. a lo<jv, <» — a — g e lu i qu' o n — a — ^voqu6 , 

l o quol — tee — iLncffiSn is g — sont — — coag — 1-es — coef f io: 



5 d 1 abaorption - dee mat e r ifcftuas, gxjsm et — de — calcul e r — g^ ux - o i 
€>t d Q l e e a ppliqu e r 5 e ye ta m e d ' g qu o tr l o n g pr6<s 6 d e nt ; . 

On notera qu'en realite le corps humain 

c 

contient trois categories pr.incipales de tissus : les 
tissus osseux, les tissue raaigres et les tissus gras, 

10 ma is qu'on n'en considdr^ que deux, en general a cause 
de difficulty £ distinguer ces. trois categories dans 
les mesures, de sorte qu'on confond volontairement les 
tissus maigres et gras. D'autres proc£d£s sont ensuite 
appliques pour distinguer leurs proportions dans les 

15 tissus raous. 

La largeur du spectre permet de disposer 
d'un nombre beaucoup plus grand <3e bandes de mesures 
qu'il ne serait n£cessaire f Y^et de les employer toutes 
pour obtenir des resultate kup p I ^menttrrres' en 

2 0 exploitant toute 1 1 information d' absorption obtenue . 
Dans 1» article "Measurement of bone mineral using a. 
multiple-energy x-ray absorptiometry tr , par J. 
Swanpalmer, R- Kullenberg, T. Hansson, Phys .Med- Biol , 
Vol 43, 1997 :pp 379-387, ou on dispose de 23 bandes de 

25 mesures et ou on corisid£re les jirois categories de 
tiSSUS y il eat propose do cotnbiner lea groupes de 
mesures trois par trois de toutes les fagons possibles 
pour obtenir 1771 (23x22x21/2x3) syst ernes de trois 
Equations a trois inconnues qui donnent autant de 

30 groupes de rasultats. Il convient alors d ! appliquer un 
critdre de choix. Les auteure conseillent de choisir 
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reliant les mesures d ' attenuation aux epaisseurs des 
materiaux. 

On notera gu'en realite le corps humain 
contient trois categories principales de tissus : les 
5 tissus osseux, les tissus maigres et les tissus gras , 
mais qu'on n'en considere que deux en general a cause 
de difficulties a distinguer ces trois categories- dans 
les mesures, de sorte qu'on confond volontairement les 
tissus maigres et gras. D'autres procedes sont ensuite 
10 appliques pour distinguer leurs proportions dans les 

. tissus mous . 

La largeur du spectre permet de disposer 
d'un nombre beaucoup plus grand de bandes de mesures 
qu'ij ne serait necessaire pour obtenir une" 7 solution, 
15 et de les employer toutes pour obtenir des resultats 
plus precis en exploitant toute 1 ' information 
d 1 absorption obtenue . Dans 1 ' article "Measurement of 
bone mineral using a multiple-energy x-ray 
absorptiometry", par J. Swanpalmer, R. Kullenberg, T. 
20 Hansson, Phys . Med . Biol , Vol 43, 1997:pp 379-387, ou on 
dispose de 23 bandes de mesures et ou on considere les 
trois categories de tissus, il est propose de combiner 
les groupes de mesures trois par trois de toutes les 
fagons possibles pour obtenir 1771 (23x22x21/2x3) 
25 systemes de trois equations a trois inconnues qui 
donnent autant de groupes de resultats. II convient 
alors d f appliquer un critere de choix. Les auteurs 
conseillent de choisir comme resultat veridique celui 
qui est a la valeur moyenne, ou a la valeur mediane, 
30 pour le parametre le plus important, qui peut etre la 
masse osseuse traversee . 



B 14155. 3/JCI 



. . MT . • 1er depot iMc^fiee : 1 0^W2 



comme r£sultat v6ridique celui qui est A la valeur 
moyenne, ou & la valeur m£diane, pour le paramdtre le 
plus important, qui peut §tre la masse osseuse 
tranrers6e - 

5 Ccmtrairement 3 ce qu ' on a conseille plus 

haut, lec auceura considerent les trois categories de 
tissus. Cela ne met pas en question la validity de leur 
methods , & condition d'expoaer le patient & une. 
intensite d ' irradiation beaucoup plus forte pour 

10 xamener les incertitudes dee meaurtB & des proportions 
acceptablee. XI serait d 1 ailleurs possible de modifier 
simplement leur methode pour I'appliquer a des rnesures 
sur deu^c tissue, ou plus g£n6ralement deux categories ; 
de raateriaux, aeulement , Toutefois, il subsiste le 

15 d€faut que cette methode ne cient pas cornpte du bruit ^ 
sur les mtsurea ec qu'elle conduit done & des resultats 
bruit 6s & e-E — i4 — apparatt quo boauoovAp — de — grottpes — 

gioouroo oont diff crcnts ; ete — cort o que lcui? r 

aombin&iflon — JV© — produit-pa g — ele- — s olutio n- — p x ecioe — aux 

2 0 cygt^taes — d'equationo oft — ±-3~& — eonfe — tst as Q Q±6o — cfc g u ' mie - 
grando — partic — de-s — calcula de ' combi n aicon — oot done 
gujette a oautioa . 

L 1 invention . a pour ob j et un 

perfect ionnement de tels proc£d£s par combinaison de 

25 reeultats nombreux et comprend une amelioration du 
critere de choix des resultats. Plus p.rdcisement , elle 
conq^me sous sa formes la plus generale un proc6d£ 
d'examen' radiologique d ' un objet ou sont considerees 
yndeux categories de materiaux, comprenant : l'emploi 

30 d'un rayonnement a spectre large ) des mesures du 
rayonnement par bandes du spectre ; des expressions 
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Contrairement a ce qu'on a conseille plus 
haut, les auteurs considerent les trois categories de 
tissus. Cela ne met pas en question la validite de leur 
methode, a condition d'exposer le patient a une 
intensite d 1 irradiation beaucoup plus forte pour 
ramener les incertitudes des mesures a des proportions 
acceptables. II serait d'ailleurs possible de modifier 
simplement leur methode pour l T appliquer a des mesures 
sur deux tissus, ou plus generalement deux categories 
de materiaux; seulement . Toutefois, il subsiste le 
t defaut que. cette methode ne tient pas., compte du bruit 
sur les mesures et qu'elle conduit done a des resultats 
bruites . 

L 1 invention a pour objet un 

perf ectionnement de tels procedes par combinaison de 
resultats nombreux et comprend une amelioration du 
critere de choix des resultats. Plus precisement, elle 
concerne sous sa forme la plus generale un procede 
d'examen radiologique d 1 un objet ou sont considerees au 
moins deux categories de materiaux, comprenant : 
l'emploi d'un rayonnement a spectre large ; des mesures 
du rayonnement par bandes du spectre ; des expressions 

{M) d'epaisseurs ou de masses des deux categories de 
materiaux traversees par le rayonnement, les 

expressions (M) etant des fonctions d'au moins deux 
des mesures (mes k ) et de coefficient (A) ; et 
1 ' application d ! un critere de choix parmi les 

expressions ( M ) pour en deduire une expression 

(M finale) jugee veridique ; caracterise en ce que le 
critere de choix comprend un combinaison (f) des 
expressions avec des facteurs de ponderation (a) , et un 
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(M) d 1 £paisseurs ou de masses des deux categories de 
materiaux travers£es par le rayormemant, les 

expressions ( M ) etanc des fonctions d'au moins deux 
des raesures (meej ec de coefficient. (A) ; et 
5 1 ' application d T uu critdre de clioix parmi le* 

expressions (M) pour en dSduire une expression 

(M finale) jugSe ' vSridique ; caractSrise en ce que le 
critere de ch.oix comprend un combinaison (f) des 
expressions avec des facteurs de pond6ration . (a) , et un 
10 calcul de factcurs de ponderation de mani£re que la A \ 

contbinaison A ait un? variation minimal d 1 aprds d e c - * <\ 

JU ^t&4^ U de^ tnesures. Cv^Xoa^ Mxr- CcxiciA&- • Wbi^ 

L * invention sera raaintenant d^crite en ' — 

r£££rence aux figures : 
15 -la figure 1 est une vue d'un spectre r 

la figure 2, ecc une repartition de 

r^sultats, 

et la figure 3 est un organigratrmi© 
r£suniant 1 1 invention, auquel on pourra se reporter dans 

20 toute la description qui suit. 

L 1 at t^nuation du rayonnernent peut €tre 
expritnee par une fonction des fipaisseurs traversSes de 
chacun de© mat6riaux d' indices x et y r ou de leurs 
masses M (densitfis par unite de surface) dans la 

25 direction du rayonnernent . Le spectre de mesureg de la 
figure l est divis€ en N bandes not£es g£n<§ralement par 
les indices i et j , Lea attenuations varieront dans 
chacujn^ des bandes en raison de coefficients 
d 1 absorption variables pour les deux materiaux. Si nous 

3 0 appelons raesi ou mesj les mesures pour une baxide 
d f energie i ou j , les masses traversees M x et My 
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calcul de facteurs de ponderation de maniere que la 
combinaison ait un bruit minimal (variance minimale en 
langage mathemat ique) calcule d'apres les bruits sur 
les mesures (variances sur les mesures) . 

L 1 invention sera maintenant decrite en 
reference aux figures : 

- la figure 1 est une vue d'un spectre, 

la figure 2 est une repartition de 

resultats , 

et la figure 3 est un organigramme 
resumant I 1 invention, auquel on pourra se reporter dans 
toute la description qui suit. 

L ' attenuation du rayonnement peut etre 
exprimee par une fonction des epaisseurs traversees de 
chacun des materiaux d' indices x et y, ou de leurs 
masses M (densites par unite de surface) dans la 
direction du rayonnement. Le spectre de mesures de la 
figure 1 est divise en N bandes notees generalement par 
les indices i et j . Les attenuations varieront dans 
chacune des bandes en raison de coefficients 
d 1 absorption variables pour les deux materiaux. Si nous 
appelons mesi ou mesj les mesures pour une bande 
d'energie i ou j, les masses traversees M x et M y 
pourront chacune etre exprimees aussi, par exemple, par 
.la formule generale 

M =Ai+A 2 . meSi+A 3 .mesj+A4 . mesf +A 5 .mes ^ +A 6 . mesi .mesj . . les 

mesures considerees dans cet exemple etant des meusures 
d 1 attenuation, on aura pour chaque canal de mesure i 
(correspondant a une bande du spectre) la relation mesi 

= In ( no y^i) °^ noi estle nombre de phatons arrivant sur 

l'objet et ni celui des photons ayant traverse I'objet. 
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pourront chacune Stre exprimees aussi, par 'example, par 
la formule generale 

r * 

M ^Ai+A 2 .meSi+Ai .mesj+A* .mes'r +A 5 .raes j +A 6 .mesi .mes-j , . Icomme 

la rion- linearity dee £onctionfl M x et M y en fonction des 
5 mesures est faible en pratique, on peut s 1 ac commode r de 

cette fonction polynomiale du deuxierne degr£ ^ qui ~ Joe. 

comprend si* coefficients A, a A 6 . 

Ces coefficients peuvent etre trouv6s dans ^ 
une 6tape d 1 6taloruiage i, travers dee etalons, appelSs 

10 pax^fois fant6mes o\i cales dans lVart f et qui consistent 
en des pieces fuiraeey d 1 epaisseurs connues et 
differentes entre ellea de materiaux simulant; par 
le\irs propri§t£s d f attenuation, les materiaux de 
l'objet qui sera ef f ect ivement a mesurer. Chacun des 

15 etalons est done soumis a \me irradiation par le 
rayonnement pencLant line longue duree qui pennet de 
require 1 r influence du bruit sur leg mesures. £e 
spectre mesur6 pour chacun des 6talons donne encore W 
mesures resultant de la decomposition du spectre en 

2 0 autant de bandes . En combinant ma in tenant deux series gjvC^- 

de mesures i et j prises pour deux bandes et pour JL^^L^J^ ' 

chacun des etalons; on recherche les coefficients A i 

pour ajuster les fonctions M x ec M y aux mesures. Dans le . 

cas present, ou il y a six coefficients A pour chacune A 

25 des deux fonctione, et ou on considers les combinaieons ^ 
d« deux bandes de mesure, lea mesures pour chacune des 
bandes devront , pour foumir une solution unique, J*^*' 
porter sur six etalons. Des etalons plus nontbreux 
pourront eiussi 6tre utilises pour amfeliorer la 

30 precision sur la solution. Une minimisation de fonction 
d'erreur sera alors appliquee. 
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Comme la non-linearite des fonctions M x et M y en 
fonction des mesures est faible en pratique, on peut 
s ' accommoder de cette fonction polynomiale du deuxieme 
degre qui comprend six coefficients Ai a A 6 . Le degre du 
5 polynome peut etre ajuste en fonction du probleme ; Par 
exemple, pour 1" analyse d'objets const itues de 
materiaux ayant des numeros atomiques plus eleves que 
celui des tissus biologiques, comme le controle non 
destructif pour I'examen d'un objet metallique . les 
10 mesures considerees dans cet exemple etant des meusures 
d f attenuation, on. aura pour chaque canal de mesure i 
(correspondant a une bande du spectre) la relation mesi 

= In ( no y .) ou noi estle nombre de phatons arrivant sur 
/ni * 

l 1 objet et ni celui des photons ayant traverse l r objet. 

15 Ces coefficients peuvent etre trouves dans 

une etape d'etalonnage a travers des etalons, appeles 
parfois fantomes ou cales dans l'art, et qui consistent 
en des pieces formees d'epaisseurs connues et 
differentes entre elles de materiaux simulant, par 

20 leurs proprietes d 1 attenuation, les materiaux * de 
1' objet qui sera ef f ect ivement a mesurer. Chacun des 
etalons est done soumis a une irradiation par le 
rayonnement pendant une longue duree qui permet de 
reduire I 1 influence du bruit sur les mesures. Le 

25 spectre mesure pour chacun des etalons donne encore N 
mesures resultant de la decomposition du spectre en 
autant de bandes . En combinant maintenant deux series 
de mesures i et j prises pour deux bandes et pour 
chacun des etalons, on recherche les coefficients A 

30 pour ajuster les fonctions M x et M y aux mesures. Dans le 
cas present, ou il y a six coefficients A pour chacune 
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Cette determination des coefficients A est 
r€p£t6a pour di verses combinaisons de pairee de 
mecurec, II en etait encore ainei dans l'article 
ant^rieur mentionn6 plus haut / cependant , *on a 
constate qu'il etait inutile d'effectuer toutes lee 

combinaisons, au nombre de — pour- exploiter 

cortipldteraent les mesures et que (N-l) combinaison 
<§taient en realite suffisantes pour recueillir toute 
1 ' information . 

Une fagon pr£f£r£e de ,p.roceder consiste a 
choisir au depart la bande de mesuree qui est la moins V-i 0 J ^ 
br-uit<Se (par exernple celle qui a le sigiial^^ 

important) et £ I'associer successivement 3 chacune; des aiA^ftH 2 ^ 
autres bandes de mesures pour les combinaisons . U "*0n^ ^ 
obtient finalement (N-l) estimations des de-ux fonctions n 
M x et M y , qu'on note M^M^.^MN -I pour chacune de cee .v^^Uc. . 
deux fonctions. 

A ce stade du proc<Sd£, les fonctions M x et 
M Yl representative* des longueurs traversees de deux 
materiaux repr^sentatif s des tissus osseux et des 
tissus mous, pourront §tre converties en £onctionSM Uy M v X 
et M w representatives des longueurs Squivalentes 
traversees des tissus osseux, des tissus maigree et des 
tissus gras en combinant lin6airement M x et M y de trois 
fa<?ons diff€rentee d£termin€es par 1 1 experience . Ce 
proc£d6 de conversion est ind£pendant de 1' invention, 
deja connu dans l'art^ et no 6 - ^ applique d l ai j kl e u r a — qu 1 & • 
l-*-4 feude d 1 ere^gt mismcs v ivaiiLa •.- • 

Comrne il n'y a pas de raison de pr<$£erer 
l'une ou 1' autre de ces estimations M , un crit£re de 
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des deux fonctions, et ou on considere les combinaisons 
de deux bandes de mesure, les mesures pour chacune des 
bandes devront, pour fournir une solution unique, 
porter sur six etalons. Des etalons plus notnbreux 
pourront aussi etre utilises pour ameliorer la 
precision sur la solution. Une minimisation de fonction 
d'erreur sera alors appliquee. 

Cette determination des coefficients A est 
repetee pour diverses combinaisons de paires de 
mesures. Il en etait encore ainsi dans 1 1 article 
anterieur mentionne plus haut ; cepehdarit, on a 
constate qu'il etait inutile d'effectuer toutes les 

Nx(N-l) 

combinaisons, au nombre de ~ , pour exploxter 

completement les mesures et que (N-l) combinaison 
etaient en realite suffisantes pour recueillir toute 
1 1 information. 

Une fagon preferee de proceder consiste a 
choisir au depart la bande de mesures qui est la moins 
bruitee (par exemple celle qui a le signal n - le 
nombre de photons a- la reception- le plus important) et 
a l'associer successivement a chacune des autres bandes 
de mesures pour les combinaisons. On obtient finalement 
(N-l) estimations des deux fonctions M Xi et M y , qu 1 on 
note M,,M 2 ,...,MV-1 pour chacune de ces deux fonctions. 

A ce stade du procede, les fonctions M x et 
M Y , representatives des longueurs traversees de deux 
materiaux representatif s des tissus osseux et des 
tissus mous, pourront etre convert ies en fonction M 0 , M v 
et M w representatives des longueurs equivalentes 
traversees des tissus osseux, des tissus maigres et des 
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choix doit Stre appliqu£ pour obtenir 1 ' estimation 

Mfinale qui sera jug6e veridique, Dana l'article 
ant^rieur, une dee expressions obtenues <§taic 
directement select ionn^e d'apr£s un critere . de 
5 classement (le resultat median) ou de moyenne das 
valeurs prises par les expressions pour un dee 

r^sulnatis . Dans 1' invention, les expressions M seront 
cotnbineeB, par example lineairement , d'apres la formule 

M finale - (ai MI) + (a2 M2)+.,+ (a N-l MN -1) tout 
10 en minimi cant la bruit ; les coefficients ai , etc. ont 

/ une somme 6gale a l 1 unite (al + a2 + .-.+ a N-l = 1), , 
Pour chaque canal de ineaure, le bruitv^suit pUokriN^o 
une loi statistigue de Poisson, dont le resultat est r 
independant pour chacun des canaux. La matrice d« 
15 covariance des N-l r<§sultats peut §tre exprimee d'apr&s 

la formule T = Y ^ ^ f 1 

La variance eur la combinalson lin£air© - 
dormant M finale s'exprime par la formule f =(a lf ^. r a N : 
i) < r . t (a x a N -i) ; cette quuantite f atteint une valeur 
2 0 optimale quand sa d6riv<§e selon toutes les variables 
est nulle, c'est-^-dire que . 1* influence du bruit est 

df/da 2 =0 ■ 



minimis^e, soit : 



La resolution numex-ique de ce syeteme 
25 fournit les coefficients al, a2, etc et M finale, 
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tissus gras en combinant lineairement M x et M y de trois 
fagons differentes determinees par 1 1 experience . Ce 
procede de conversion est independant de 1' invention, 
de j a connu dans 1 1 art , 
5 Comme il n'y a pas de raison de preferer 

1 ' une ou 1 1 autre de ces estimations M , un critere de 
choix doit etre applique pour obtenir l 1 estimation 

M finale qui sera jugee veridique . Dans I 1 article 
anterieur, une des expressions obtenues etait 
10 directement . selectionnee d'apres un critere de 
- classement (le resultat median) 4 bu de moyenne des 
valeurs prises par les expressions pour un des 

resultats. Dans 1* invention, les expressions M seront 
combinees, par exemple lineairement, d'apres la formule 

15 M finale = (al Afl) + (a2 M 2 ) +...+ (a N-l M N -1). tout 
en minimisant le bruit ; les coefficients a X/ etc. ont 
une somme egale a 1 ' unite (al + a2 + a N- 1 .= 1). 

Pour chaque canal de mesure, le bruit susr 
le nombre de. photons suit une loi statist ique • de 

20 Poisson, dont le resultat est independant pour chacun 
des canaux. La mat rice de covariance des N-l resultats 
peut etre exprimee d'apres la formule 

ij dmes k dmes k N k 

La variance sur la combinaison lineaire 
25 donnant M finale s 1 exprime par la formule f =(ai,...,a N _ 
i) -T. fc (a lf a N -i) ; cette quantite f atteint une valeur 
optimale quand sa derivee selon toutes les variables 
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c'est-a-dire les masses traversees des deux categories 
de materiaux, 

Ce soat ces operations qui sont accomplies 
dans 1 ' invention ; il est important de* rcmarquer que 
5 les mesures (mask) employees dans la formule dormant T y 
et exieuite sent les mesures effectuees 3. travers 
I'objet & £tudier lui-meme, mais pas cellee qui ont etS 
faiths a 1'^talonnage pour determiner les coefficients 

A des- fonctiona M „ 

10 Li 1 invention permet d'£chapper au dilemme de 

proced£s anterieurs dans lesquelsXon renon<?ait & une 
partie des Energies de mesure ou au contraire on 
acccptaic des mesures peu distinctes : elle utilise en 
effet tout le spectre, mais lo divieant en bandes aseez 

15 nombreuses pour que la mesure de chacune puisse etre 
comparee de fa<?on utile a d'autres mesures, faites «ur 
des bandes eloign€es . Elle' permet ainsi, entre autr^s 
avantages, d'etudier aussi bien des organismes maigres 
que gras , 

20 Enfin, 1 ' invention peut etre g£n£ralis£e £ 

un nombre de mat^riaii/plus grand que deirx, ce qui peut X 
etre interessant notammenc pour les precedes d'imagerie 
a produit de contraste, ou troia variables doi.vent etre 
considergea ) ec elle peut^ etre appliquee en 

25 considerant des combinaisoris autres que lin6aires des 

fonctions M , &j^- *x<amp lm < w appliqu a t^ — l e proo < 5d6 A CP 
• (toalyoo — — aompooant — Principal) — conna dana — 1 r art — efe — 
n ^ntionn6 — p&z — set ample — daas — 5 b ^articl e — de — M Aitr ejaan, 

Per ion Ouiideymjiuii e^fc im;icul6 ilHon i-m ^t- i w 

3 0 chfami cal identifi ca tion — u s in g — as — X ray — trah e miogi e n 

function obtain e d with th o multi -e nergy ro g ttod" , — grcNg , 
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est nulle, c'est-a-dire que 1 ' influence du bruit est 

'9/73(2, =0 
df/da 2 =0 



minimisee, soit : < 



df/da N ^=Q 



La resolution numerique de ce systeme 

5 fournit les coefficients al, a2 , etc et M finale, 
c'est-a-dire les masses traversees des deux categories 
de materiaux. 

. Ce sont ces operations qui sont accomplies 
dans l r invention ; il est important de remarquer que 
10. les mesures (mes k ) employees dans la formule donnant T tJ 

et ensuite sont les mesures effectuees a travers 
l ! objet a etudier lui-meme, mais pas celles qui ont ete 
faites a I'etalonnage pour determiner les coefficients 

A des fonctions M . 

15 L 1 invention permet d'echapper au dilemme de 

procedes anterieurs dans lesquels on renon^ait a une 
partie des energies de mesure ou au contraire on 
acceptait des mesures peu distinctes : elle utilise en 
effet tout le spectre, mais en le divisant en bandes 

2 0 assez norabreuses pour que la mesure de chacune puisse 
etre comparee de fagon utile a d'autres mesures, faites 
sur des bandes eloignees. Elle permet ainsi, entre 
autres avantages, d ! etudier aussi bien. des organismes 
maigres que gras . 

25 Enfin, 1 1 invention peut etre generalisee a 

un nombre de materiaux plus grand que deux, ce qui peut 
etre interessant notamment pour les procedes d f imagerie 
a produit de contraste, ou trois variables doivent etre 
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S PIE ^ol - — 3446 ; p , 134 -15 2 , Hard .X-r a y — and g a mm a -r ay 
do tr o otror physios a nd - a p p lication s , — 6dit§ par F , — Patrick 
Poty Qt Richard B . Hoove r, 
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considerees ou pour controle des bagages (recherche 
d'explosifs) ; et elle peut etre appliquee en 
considerant des combinaisons autres que lineaires des 

f one ti oris . M . . . , . 

Des application de 1 1 invention sont : 
. 1 1 osteodensitometrie 

- obtention de la densite de masse 

osseuse 

obtention de la composition 
corporelle (masse grasse, masse maigre) 

. le controle agro-alimentaire, par exemple 
la detection de fragments d'os dans de la viande ou la 
detection de morceaux de verres dans des plats 
cuisines , 

le controle des bagages : recherche 
d T explosifs, de produits illicites (armes, aliments, 
drogue...) 
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REVINDICATIONS 



10 



15 



20 



25 



1) Precede d'e^amen radiologique d'un objet 
oii sont considerees V deux categories de materiaux, 



a spectre 
bandes du 



comprenant : l'emploi d'un rayonneraent 
large ; des mesures du rayonnement par 

spectre ; des expressions (M) d'epaisseurs ou de 
nuflces des deux categories de materiaux travers6cs par 

le rayonnement, les expressions (M) £tant des 
fonctions d 1 au moins d^ux des mesures (mes k ) et de 
coefficients (A) ; et 1 T - application d'un crit§re de 
ch-oix parmi les expressions ( M ) pour 6n deduire une 

expression {M finale) jug6e v€ridique ? caract£ris£ en 

ce que le critere de choix comprend une combinaison (f) 

dee expressions avec des facteura de pond^ration (a) , 

et un calcul des facteurs de pond€ration de maniere que. . 

ajuiv SrujuJb GcUculeT ' 

la combinaison ait \m® — variaclon minimal* d'apr^s de«* Je. 
^xut He/* ' 

2) Proc6de selon la. revendication 1, 
caract^rise en ce que la combinaison (f) des 
expressions eat lineaire. 

3) Proced6 selon 
caracteris^ en ce que la variation de la combinaison 



i-dfe revendicationS 2, 



est calcul£e avec une matrice de covariance des 

expressions ( M ) . 

4) Procede selon l'une quelconque des 

revendications 1 & 3 , caract6rise en ce que les 

expressions (M ) sont en nombre <§gal a celui des 

bandes, moins un, et etablies coujours avec une des 
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RE VEND I CAT I ONS 

1) Procede d'examen radiologique d ! un objet 
ou sont considerees au moins deux categories de 
materiaux, comprenant : 1 1 emploi d'un rayonnement a 
spectre large ; des mesures du rayonnement par bandes 

du spectre ; des expressions (M) d { epaisseurs ou de 
masses des deux categories de materiaux traversees par 

le rayonnement, , les expressions ( M ) etant des 
fonctions d'au moins deux des mesures (mes k ) et de 
coefficients (A) ; et 1 ' application * d f un critere de 

choix parmi les expressions (M) pour en deduire une 
expression (M finale) jugee veridique ; caracterise en 
ce que le critere de choix comprend une combinaison (f) 
des expressions avec des facteurs de ponderation (a), 
et un calcul des facteurs de ponderation de maniere que 
la combinaison ait un bruit minimal calcule d'apres un 
bruit sur les mesures. 

2) Procede selon la revendication 1, 
caracterise en ce que la combinaison (f) des 
expressions est lineaire. 

3) Procede selon 1 ' une quelconque des 
revendications 1 ou 2, caracterise en ce que la 
variation de la combinaison est calculee avec une 

matrice de covariance ( V i; ) des expressions ( M ) . 

4) Procede selon Pune quelconque des 
revendications 1 a 3, caracterise en ce que les 

expressions ( M ) sont en nombre egal a celui des 
bandes, moins un, et etablies toujours avec une des 
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bandes (io) et r reapectiveitient, chacune des autree 
bandes , 

5) Proc^de selon 1'une quelconqua des 
revendi cat ions 1 & 4, caract;£ris£ en ce que; lea 
coefficients (A) sont determines a une §tape pr£alabla 
de calibration. 

<^V^ ^ ^^^^ 
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bandes (io) et, respect ivement , chacune des autres 
bandes . 

5) Procede selon 1'une quelconque des 

revendications 1 a 4, caracterise en ce que les 

- * - . . 

coefficients (A) sont determines a une etape prealable 
de calibration. 

6) Procede selon l'une quelconque des 
revendications precedsente, caracterise en ce qu'il est 
applique a 1 ' osteodensitometrie . 

7) Procede selon 1 1 une quelconque des 
revendications 1 a 5, caracterise en ce qu'il est 
applique a des controles agro-alimentaires . 

8) Procede selon l'une quelconque des 
revendications 1 a 5, caracterise en ce qu'il est 
applique au controle de bagages . 
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MESURES A TRAVERS DES ETALONS 
POUR UN SPECTRE A N BANDES (ij) 



CALCUL DES COEFFICIENTS (A) 
DES FONCTIONS DONNANT LES EPAISSEURS 
TRAVERSEES (CONNUES) D'APRES 
LES MESURES (mes- ( .mesj ) 



MESURES A TRAVERS L'OBJET A ETUDIER 



j ESTIMATIONS (M) DES EPAISSEURS 
TRAVERSEES EN APPLIQUANT 
LES FONCTIONS D'ETALONNAGE 



COMBINAISON* (f) DES ESTIMATIONS 
AVEC DES FACTEURS DE PONDERATION (a) 



AJUSTEMENT DES FACTEURS DE 
PONDERATION POUR MINIMISER 
UNE VARIATION DE COMBINAISON 
EN FONCTION DES MESURES 



ESTIMATION FINALE DES EPAISSEURS 
EN APPLIQUANT LA COMBINAISON (f) <- 
ET LE^FACTEURJAJUSTES 



FIG. 3 
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MESURES A TRAVERS DES ETALONS 
POUR UN SPECTRE A N BANDES (i,j) 



CALCUL DES COEFFICIENTS (A) 
DES FONCTIONS DONNANT LES EPAISSEURS 
TRAVERSEES (CONNUES) D'APRES 
LES MESURES (mesj ,mesj ) 



MESURES A TRAVERS L'OBJET A ETUDIER 



ESTIMATIONS (M) DES EPAISSEURS 
TRAVERSEES EN APPLIQUANT 
LES FONCTIONS D'ETALONNAGE 



COMBINAISON (f) DES ESTIMATIONS 
AVEC DES FACTEURS DE PONDERATION (a) 



AJUSTEMENT DES FACTEURS DE 
PONDERATION POUR MINIMISER 
UNE VARIATION DE COMBINAISON 
EN FONCTION DES MESURES 



ESTIMATION FINALE DES EPAISSEURS 
EN APPLIQUANT LA COMBINAISON (f) 
ET LES FACTEURS AJUSTES 



FIG. 3 
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